Two-photon absorption in potassium niobate 
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We report measurements of thermal self-locking of a Fabry-Perot cavity containing a potassium 
niobate (KNbOa) crystal. We develop a method to determine linear and nonlinear optical absorp- 
tion coefficients in intracavity crystals by detailed analysis of the transmission lineshapes. These 
lineshapes are typical of optical bistability in thermally loaded cavities. For our crystal, we deter- 
mine the one-photon absorption coefficient at 846 nm to be a = (0.0034 ± 0.0022) m _1 and the 
two-photon absorption coefficient at 846 nm to be /3 = (3.2 ± 0.5) x 10 -11 m/W and the one-photon 
absorption coefficient at 423 nm to be (f 3 ± 2) m -1 . We also address the issue of blue-light- induced- 
infrared- absorption (BLIIRA), and determine a coefficient for this excited state absorption process. 
Our method is particularly well suited to bulk absorption measurements where absorption is small 
compared to scattering. We also report new measurements of the temperature dependence of the 
index of refraction at 846 nm, and compare to values in the literature. 



I. INTRODUCTION 

For cw laser second harmonic generation in a build-up 
cavity, thermal issues in the nonlinear crystal inevitably 
come into play. This is particularly true at high laser 
powers. In typical applications, a Fabry-Perot build-up 
cavity increases the fundamental laser intensity by one 
or two orders of magnitude compared to the single-pass 
intensity. The nonlinear crystal absorbs a fraction of 
the fundamental power, heats up, and changes the cavity 
properties. Thermal lensing changes the cavity's spatial 
mode. Nonlinear absorption changes the cavity's finesse. 
Thermal expansion and changes in the crystal's index of 
refraction change the cavity's optical path length. 

These effects are well known, and have been observed 
by several groups working in this field (for example, 
|], g]). Thermally- induced changes in the build-up cav- 
ity sometimes limit the second-harmonic generation ef- 
ficiency In many cases, they are nuisances that can be 
worked around. However, these changes can be used as 
a sensitive diagnostic tool for accurately characterizing 
nonlinear crystal properties. 

For gases and atomic vapors, measuring optical prop- 
erties of materials inside Fabry-Perot cavities has a long 
history particularly in relation to optical bistability, 
saturable absorption, and index of refraction measure- 
ments. Thermal and collisional properties of gases have 
also been measured in this way Q, as well as optical 
properties of cavity mirrors ||. However, measurements 
of optical properties of crystals in Fabry-Perot cavities is 
relatively new ||. 

In some ways, Fabry-Perot measurements are simi- 
lar to laser calorimetry 0. In those calorimetry mea- 
surements, a laser passes through a thermally isolated 
absorbing crystal and heats it slightly. From the tem- 
perature rise, laser power, thermal mass, and thermal 
conductivity it is possible to determine an optical ab- 
sorption coefficient for the crystal. Comparing calorime- 
try to our Fabry-Perot measurements, both are pho- 
tothermal methods. Both can have high sensitivities to 
small absorptions. In contrast to calorimetry, our "ther- 



mometer" is the power-induced shift in the cavity res- 
onance frequency. Fabry-Perot measurements can typi- 
cally reach much higher intensities, making nonlinear ef- 
fects more apparent. Furthermore, frequency-based tech- 
niques, such as the one presented in this paper, usually 
offer a higher level of precision. 

This paper presents a detailed study of the thermal 
response of KNb03 inside a Fabry-Perot build-up cav- 
ity. We use this response to determine a new value for 
the optical absorption coefficients in KNb03 at 846 nm 
and 423 nm, and to measure the effects of bluc-light- 
induced infrared absorption (BLIIRA). We also measure 
the temperature-dependence of the index of refraction for 
light polarized along the b and c crystal axes. These mea- 
surements are compared to values from the literature. 



II. EXPERIMENTAL SETUP 

Our experimental setup is typical for cw second- 
harmonic generation in an external build-up cavity (see 
Figure [l]). This laser system is designed for a calcium 
laser cooling experiment at 423 nm. It produces typi- 
cally 75 mW at 423 nm, measured outside the cavity. 

The laser system consists of a single frequency ex- 
tended cavity diode laser at 846 nm amplified in a single- 
pass through a tapered laser diode || . The extended cav- 
ity diode laser is frequency-stabilized relative to a sepa- 
rate optical cavity using the Pound-Drever method || to 
provide both short and long-term stability. After optical 
isolation, polarization correction, and mode matching, we 
inject 200 mW into the Gaussian mode of a four-mirror 
folded (bowtie) Fabry-Perot frequency doubling cavity. 
Higher-order modes in the cavity are less than 5% of the 
Gaussian mode. The input coupler reflectivity is 96.6%, 
and the other three mirror reflectivities are > 99.5% in 
the infrared. The radius of curvature for the two curved 
mirrors is 100 mm, and they are separated by a distance 
of 117 mm. The angle of incidence on the curved mirrors 
is 6.5°, minimizing the cllipticity of the beam waist inside 
the crystal (41.1 and 40.7 ^m in the tangential and saggi- 
tal planes, as determined from the cavity geometry). The 



2 



ECDL 




Amplifier 






FIG. 1: Schematic diagram of the experimental layout. ECDL 
= extended cavity diode laser, I = optical isolator, A/2 = half 
wave plate, A/4 = quarter wave plate, CI = frequency dou- 
bling "bowtie" cavity, C2 = stable optical cavity, D = photo- 
diode detectors, PBS = polarizing beam-splitter cube, BS1 = 
uncoated quartz window, BS2 = dichroic beam splitter. The 
short heavy lines are mirrors. 



round-trip cavity length is 618 mm. The KNbC>3 crys- 
tal jlO| is a-cut and antireflection-coated (R <0.25%) at 
846 nm. The non-critical phase-matching temperature is 
— 11.5°C, for light polarized along the 6-axis. The crystal 
is mounted on a Peltier-cooled copper block. An AD590 
temperature sensor and a thermister embedded in the 
copper mount monitor the crystal temperature. The en- 
tire bowtie cavity is enclosed in a sealed aluminum box, 
filled with dry oxygen at atmospheric pressure. For low 
laser powers, and when the crystal temperature is far 
from the phase matching temperature, the cavity finesse 
with the crystal in place is 124. Without the crystal, the 
cavity finesse is 140. 



III. TRANSMISSION LINESHAPE WITHOUT 
BLUE LIGHT 

When laser light incident on the bowtie cavity is polar- 
ized along the c-axis, it is impossible to meet the phase- 
matching condition for non-critical second harmonic gen- 
eration at A = 846 nm, and no blue light is generated. 
Previous measurements in KNb03 have shown that even 
small amounts of blue light dramatically change the in- 
frared absorption properties jll], |lj, |l3) . We address this 
blue-light-induced-infrared-absorption (BLIIRA) later in 
the paper. However, in this section we discuss infrared 
absorption in the absence of blue light, i.e.: when the 
incident laser beam is polarized along the c-axis. 

One of the flat mirrors in the bowtie cavity is mounted 
on a piezoelectric crystal, and we use this crystal to 
change the optical length of the cavity. A small amount 
of infrared light is transmitted through the cavity mir- 
rors, and we measure this light while sweeping the 



cavity length. When the piezo-mounted mirror passes 
through the position for cavity resonance (i.e.: when the 
roundtrip pathlength in the cavity is an integer multi- 
ple of the laser wavelength), the power circulating in the 
cavity, and transmitted through the cavity end mirror, 
increases. The transmission lineshape is Lorenzian in the 
absence of disturbing effects: 



P{v) 



Po 



1 + F (y - Vq) 



(1) 



where i>o is the resonant frequency of the cavity, 2/y/F is 
the full width at half-maximum (FWHM), and P is the 
infra-red power at the peak of the Lorenzian. 

The KNb03 crystal inside the bowtie cavity absorbs 
a small amount of the laser power and heats up. This 
changes the crystal's length and index of refraction, 
which in turn changes the optical path length of the 
bowtie cavity. We model the overall effect by replacing 
vq with uq + tjP/Po. Accordingly, Equation |l| becomes 



V = 



i + f (v — m/Y 



(2) 



where y = P/Pq is the normalized transmission, and 
the frequency v is measured relative to the unshifted 
frequency vq. Transmission lineshapes of this kind are 
typical of optically bistable devices. A similar treatment 
of this effect is given in the references . 

The peak in the transmission curve occurs when y = 1 
and v — -q. Figure ^| shows the transmission lineshape for 
both low and high laser powers. On the same figure we 
plot a fit of the data to Equation^. The asymmetry in the 
transmission lineshape at high powers is due to thermal 
changes in the optical pathlength. When sweeping the 
cavity from low to high frequencies as in Figure |^, the 
power-dependent cavity resonance is effectively "dragged 
along" from v — v to v = uuq + rj as the cavity length 
changes. This is known as "thermal self-locking" Qj. The 
value of r\ is called the "self-locking range," and it is 
related to the "thermo-optic constant" of reference || . 

For a given power incident on the cavity, we determine 
the self-locking range r\ using a nonlinear least-squares 
fit of the measured transmission lineshape to Equation 
||. For the highest laser intensities, it appears that a 
better overall fit is obtained by replacing r\y in Equa- 
tion ^| by a higher order polynomial. However, the fit 
in these high intensity data is hampered because the fit 
function is not single-valued, as can be seen in Figure 
||. For these measurements, we determine the values of 
•q visually by finding the frequency at the peak of the 
transmission curves. 

Our values for 77 as a function of infrared power in 
the bowtie cavity are shown in Figure ^a. The error 
bars indicate the 1 a statistical spread in 5 measurements 
at each power setting. For linear absorption, the cavity 
shift, 77, should be a linear function of the circulating 
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FIG. 2: Infrared transmission lineshape for the bowtie cavity 
in the absence of blue light. The points are the measurements. 
The solid line is a fit to the data (see Equation [|). The 
infrared light is polarized along the c-axis. (a) Low power, (b) 
High power. For high powers, the fit function is not single- 
valued above the cavity resonant frequency. 
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FIG. 3: (a) Self-locking range, rj, versus infrared detector 
signal, (b) The data from j3|a with temperature and power 
calibrations. The solid line is a weighted fit to the data, T = 
ciP + c 2 P 2 . 



ity shifts at a rate of 231 MHz °C~ 1 for light polarized 
along the 6-axis and 331 MHz °C~ 1 for light polarized 
along the c-axis. 

We calibrate our photodiode signal using a power me- 
ter, and measure the efficiency of the optical path to con- 
vert this IR signal to a power circulating inside the cavity. 
The conversion from cavity shift to temperature rise can 
also be derived from measurements of the change in the 
index of refraction with temperature. From the published 
values of dn/dT |Q and the geometry of our bowtie cav- 
ity, we calculated what the temperature-dependent shift 
ought to be. The calculation agrees with our data for the 
6-axis. However, we find a small difference between the 
calculation and our data for the c-axis, as discussed later 
in this paper. 

Using our measurements, we convert 77 into AT, the 
difference in temperature between the crystal axis where 
the laser propagates and the crystal wall which is kept 
constant by our temperature control circuit. Figure |3p 
shows a plot of the temperature rise as a function of the 
infrared power circulating inside the cavity. It is rela- 
tively simple to convert this temperature rise into a mea- 
sure of the optical absorption coefficient of KNbC>3 by 
solving the heat transfer equation, which we do below. 

IV. SOLVING THE HEAT TRANSFER 
EQUATION 

We solve the heat transfer equation in cylindrical co- 
ordinates. The laser beam is in the TEMoo Gaussian 
mode, with a 41 /im waist. The confocal parameter of 
the beam is 6.2 mm, and at the face of the crystal, the 
beam waist is only 7% larger than in the center of the 
crystal. So we approximate the laser beam as axially 
symmetric and Gaussian with a 41 /im waist along the 
entire 10 mm length of the crystal. Our crystal has a 
square cross section, 3 mm on a side. However, because 
it is so much larger than the laser beam waist, we as- 
sume that the crystal is cylindrical also, with a 1.5 mm 
radius. With these assumptions, the equilibrium temper- 
ature inside the crystal is determined by solving the heat 
equation in the form: 



power. However our measurements indicate a quadratic 
dependence of the cavity shift on the circulating power, 
suggesting that the dominant absorption mechanism is 
two-photon absorption in the infrared. 

We convert 77 into a temperature rise inside the crys- 
tal by measuring the temperature change necessary to 
shift the cavity resonance frequency by one free spectral 
range. In these calibration measurements, we attenuate 
the incident laser power to 10 /iW to avoid laser-heating 
in the crystal. The cavity mirror positions remain fixed, 
and the laser frequency is constant. Using the Peltier- 
device, we change the temperature of the crystal by sev- 
eral degrees while monitoring the cavity transmission. At 
T = — 11.5°C, the resonant frequency of the bowtie cav- 



r dr \ dr ) K ' 

where u(r) is the temperature as a function of radius, G 
is the power density, and K is the thermal conductivity. 
Following Sutherland the absorbed power density 
inside the crystal can be written as 

G = al + f3I 2 (4) 

where I = (2P/-Ka 2 ) exp(— 2r 2 /a 2 ) is the laser intensity, 
P is the laser power incident on the crystal, a is the 
Gaussian beam waist, a is the linear optical absorption 
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coefficient, and j3 is the two-photon absorption (TPA) 
coefficient. 

We can simplify Equation by letting r = ay/^/2, 
A = (Pa/2irK), and B = (P^(3/2Tr 2 a 2 K). With these 
substitutions, Equation || integrates once to 



p = - [oM-y 2 ) - 1] + - [cx P (-2 y 2 ) - I] , (5) 
ay y y 

where the constant of integration is chosen to avoid a 
singularity as y — > 0. The solution to this differential 
equation is 



u{y) 



A 



Ei (y 2 ) - ^ ln (^ )~2 El (2y2) ~ 2 Hy2) + ° 

_ (6) 

where E n (x) is the exponential integral flqj , and C is 
the constant of integration. This constant is determined 
by requiring the temperature at the wall of the crystal 
to be T c (i.e.: u(r = b) = T c ), which we control in the 
experiment. With this constraint, we can evaluate the 
temperature rise on axis inside the crystal. The series 
expansion for the exponential integral has a logarithmic 
term that exactly cancels the ln(y 2 ) terms in Equation 0. 



AT = 



u(0)-T c = ^ 

a I 
4nK n I a 2 



In 
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a- 



2b 2 



+ 7 



(7) 



where 7 = .57721 ... is Euler's constant. 

We can analyze the data in Figure || to extract the lin- 
ear and two-photon absorption coefficients. A weighted 
fit to the function T(P) = c\P + C2P 2 , gives c\ — 
(5.7±3.6)xlfr 4o C W" 1 andc 2 = (1.09±0.12) x 1CT 3 °C 
W~ 2 . The uncertainties in these numbers are the statisti- 
cal uncertainty from the weighted fit to the data. Taking 
the thermal conductivity from H K = 4.0 W m _1 °C~ 1 , 
a = 41 x 1CT 6 m, b = 1.5 x KT^m, and 7 = 0.57721. We 
determine the one-photon (linear) absorption coefficient 
to be a = (0.0034 ± 0.0021) m" 1 and the two-photon 
absorption coefficient to be j3 = (3.16 ± 0.35) x 10~ n 
m W _1 . The uncertainties in these absorption coeffi- 
cients are only statistical. We have a systematic uncer- 
tainty in the power measurements of 10%. Our assump- 
tion of cylindrical symmetry in the crystal for solving 
the heat equation also introduces an error. The error 
is known, but related to this is the fact that our laser 
does not propagate exactly down the center of the crys- 
tal. These errors show up in the ln(& 2 ) term in Equa- 
tion 0, and probably add ~ 8% uncertainty to the mea- 
surements. Our best numbers for the absorption coef- 
ficients are therefore a = (0.0034 ± 0.0022) m^ 1 and 
(3 = (3.2 ± 0.5) x lO" 11 m/W. These numbers do not re- 
flect the uncertainty in the thermal conductivity, which 
we do not know. To our knowledge, this is the first de- 
termination of the two-photon absorption coefficient at 
any wavelength in KNb03. 

Earlier work on the linear absorption coefficient for 
KNb0 3 § found a = 0.001cm- 1 = 0.1m- 1 at 860 ran. 
This value is much higher than this work. It may be that 
some of the absorption measured in previous work was 
two-photon absorption. That work used laser calorimetry 
to measure the absorption coefficient at specific laser lines 



from 457 nm to 1064 nm. At 860 nm, they focused a 900 
mW laser beam into a 5 mm crystal, and it is likely that 
the intensity was high enough for multiphoton absorption 
to be important. 

Our method for determining the one and two-photon 
absorption coefficients is quite general. It is independent 
of measurements of light scattering inside the crystal. 
It is also independent of reflection measurements at the 
crystal faces. This method is particularly well suited to 
measurements of optical absorption where the absorption 
coefficients are comparable to or smaller than scattering 
and reflection coefficients. 



V. TRANSMISSION LINESHAPE WITH BLUE 
LIGHT 

When infrared light incident on the bowtie cavity is 
polarized along the 6-axis of the KNb03 crystal, it is 
possible to meet the phase-matching conditions required 
for generating the second harmonic at 423 nm. The pres- 
ence of blue light in the crystal dramatically alters the 
transmission lineshape. Not only is the blue light itself 
absorbed by the crystal, it also significantly increases the 
absorption of the infrared light. This is called "blue- light- 
induced-infrared-absorption" (BLIIRA), and it has been 
studied at length in the literature 0, jll|, [l2|, Pre- 
vious work has found BLIIRA to be significant at blue 
light intensities down to 7 x 10~ 4 W/cm 2 . BLIIRA is 
minimized at longer wavelengths and at higher crystal 
temperatures [ p^[ . However, for our calcium work, we 
require high powers at 423 nm. 

As before, we fix the laser frequency and scan the cav- 
ity length through the resonance condition. Both blue 
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FIG. 4: (a) Cavity transmission lineshape in the infrared at 
846 nm when generating the second harmonic at 423 nm. This 
signal is taken from the highest intensities in the experiment. 
Note that for low powers, the FWHM is 5 MHz. (b) The sec- 
ond harmonic signal for the same settings as ^a. The shoulder 
to the left of the signal is a Maker fringe (see text), (c) Peak 
second harmonic power as a function of peak infrared power. 
For all but the highest two power measurements, the peak 
blue power is given by Puue = 36.2kW _1 Pf r . 



light at 423 nm and a small amount of infrared light at 
846 nm exits the cavity. We separate these wavelengths 
using dichroic mirrors, which transmit 95% in the blue 
and reflect 99.5% in the infrared. To make sure that no 
blue light reaches the infrared laser beam detector, we 
use four of these mirrors in series in the infrared beam 
path after the doubling cavity (see Figure 0). 

In these experiments, we do not independently control 
the infrared and blue beam intensities. Rather, we opti- 
mize the crystal temperature to maximize the blue light 
production for each infrared power setting in the steady 
state. Losses inside the doubling cavity limit the maxi- 
mum infrared power circulating inside and therefore the 
maximum blue light produced. 

The transmission lineshape for the cavity when blue 
light is present is shown in Figure |]a for our highest IR 
powers. The temperature rise inside the crystal is signif- 
icant. At maximum power, as the cavity approaches the 
resonance condition, only a small amount of blue light 
is initially generated because the crystal temperature is 
too low to meet the phase matching condition for sec- 
ond harmonic generation exactly (see Figure ||b) . Closer 
to the cavity resonance condition, the circulating power 
increases, the crystal heats up, and more blue light is gen- 
erated, which heats the crystal even more. This positive 
feedback continues, and the cavity demonstrates thermal 
self- locking for up to 90 MHz, as shown in Figure ^a. 

The blue light signal in Figure |i|b has a shoulder on 
the low- frequency side of the maximum. This is a Maker 
fringe. As stated previously, we initially optimize the 
crystal temperature for optimum blue light production 
in the steady state. At the highest intensities, the tem- 
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FIG. 5: (a) Self-locking range, r], in the presence of blue light 
versus infrared detector signal, excluding the highest power 
data where the SHG process appears to be saturating, (b) 
The data from [E^, with the temperature and power calibra- 
tions, and with the temperature rise from the infrared light 
subtracted off. The solid line is a weighted fit to the data (see 
text). 

perature shift (~ 0.4°C) due to thermal self-locking of 
the cavity is a few times the temperature phase match- 
ing bandwidth. 

Because of these complicated thermal conditions, and 
especially because the blue light production is not con- 
stant across the transmission lineshape, we are reluctant 
to fit the lineshape to any model for the highest laser 
intensities. Instead, we find the self-locking range graph- 
ically. It is the falling step at the far right hand side of 
the infrared transmission peak shown in Figure |. This 
choice is valid as long as there is no significant thermally- 
induced change in the cavity coupling efficiency. We can 
monitor these changes by measuring the light reflected 
from the cavity. For the relatively modest intensities 
in this study, the minimum reflected light (which corre- 
sponds to the maximum light inside the cavity) changes 
by only a few percent, and only at the highest inten- 
sity measurements. The statistical scatter in the data 
is larger than this, and we neglect this small systematic 
error. 

At lower laser intensities, where BLIIRA and other ab- 
sorption processes are less severe, we use a least-squares 
method similar to Section 3 to find the self-locking range 
rj. The self-locking range is plotted as a function of in- 
frared laser power in Figure ||a. Note that except for the 
two highest power measurements, the frequency doubling 
process is not saturated (see Figure |]c) . 

With blue light present in the crystal, the temperature 
rise inside the crystal has three sources: infrared absorp- 
tion (one and two photon), blue light absorption, and 
BLIIRA. We can easily subtract off the infrared contri- 
bution to the self-locking range using the data in Figure 
|as a look-up-table. The remaining contribution to the 
self-locking range is due to temperature rise from blue 
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light absorption and BLIIRA. 

Because the blue light intensity and the temperature 
are not constant along the length of the axis of the crys- 
tal, we must be careful about how we define the temper- 
ature rise and the blue light intensity. We will define the 
absorption coefficient at 423 nm, S, by the equation 



AP 



hSdV = 5P b L 



(8) 



where lb is the blue light intensity inside the crystal, 
Pi, is the measured blue power, and L is the length of 
the crystal. In the non-depleted pump plane-wave ap- 
proximation, the blue light intensity grows quadratically 
with distance z in the crystal. With the definition from 
Equation ||, the blue light intensity can be written as 
Ib(r,z) — (12P b z 2 / ira 2 L 2 ) exp(— 4r 2 /a 2 ), where a is the 
Gaussian waist for the infrared beam in the crystal and 
we explicitly assume that the Gaussian waist for the blue 
beam is a/^/2. 

The heat equation for blue light absorption and BLI- 
IRA is written as 



1 d ( du 
r dr \ dr 

= -z 2 (Ce- Ar ^ aA +De- 6r ' Z / a ' 



o u 

dz 2 



(9) 



The first term on the right hand side of the equation is the 
blue light absorption term where C = 12Pf,6 / '(tto 2 L 2 K) . 
The second term on the right hand side of the equation is 
the BLIIRA term, where D = 24P 6 P„£ / (tt 2 a 4 L 2 K) , P lr 
is the measured infrared power, assumed to be constant 
along the crystal length, and £ is the BLIIRA absorption 
coefficient. As before, u(r, z) is the temperature inside 
the crystal and K is the thermal conductivity. 
The solution of this equation for r = is 
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35 
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-P b P ir , (10) 



but we have to be careful about how we define the tem- 
perature rise. In Equation |^, the radial derivatives are 
much larger than the longitudinal derivatives, and the 
d 2 u/dz 2 term can be treated perturbatively. This ap- 
proximation produces a temperature function that in- 
creases quadratically with z in the crystal. 

We can write the temperature rise on axis as AT(z) = 
T m z 2 /L 2 , where T m is the maximum temperature rise 
in the crystal. We can calculate the change in the opti- 
cal pathlength inside the crystal with this temperature 
profile. It is equal to the change in pathlength due to an 
equivalent temperature AT' = T m /3. More expressly, 
we can write T{z) = T m z 2 /L 2 = 3T'z 2 /L 2 , and rewrite 
Equation ITol as 



r 



4nK 
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(11) 



Knowing this, we can convert the measured self-locking 
range of the cavity in the presence of blue light into a 
temperature rise. The data is plotted in Figure |5|. For 
this data, the "temperature rise" is the average tempera- 
ture rise in the crystal, i.e.: T". The "blue power" is the 
measured blue power, i.e.: the power at the end of the 
crystal, Pb- 

Again, we fit a function of the form 



AT = ai P b + a 2 P b P„ 



(12) 



to the data in Figure |5|b using a multivariable least- 
squares method and find ci\ = 2.41°C W _1 and a 2 — 
1.04°C W~ 2 . Using equation 11, we determine the 
absorption coefficients 5 — 13.3 m _1 and £ = 2.2 x 
10 _8 m/W. Apparently, there are no measurements in 
the literature of the KNbO*3 absorption coefficient at 423 
nm. However, there are measurements close to this wave- 
length [Q, and a value of 13.3% per cm is reasonable. For 
this determination of the absorption coefficient, we have 
the same kinds of uncertainties as before, namely in the 
beam position in the crystal and in the power measure- 
ment. We estimate the error in this measurement to be 
15 %. Our final number for the absorption coefficient at 
423 nm is S = (13 ±2) m" 1 . 

This uncertainty estimate is perhaps conservative. We 
have a second KNb03 crystal. In a simple transmission 
measurement at 423 nm, using the blue light after the 
frequency doubling cavity, we find the linear absorption 
coefficient to be 13% per cm, in agreement with the work 
above. 

We verified the validity of our perturbative approxima- 
tion by solving Equation |^ numerically. In the numerical 
solution, the temperature profile is quadratic in z ex- 
cept at the very end of the crystal, where is flattens out 
slightly to meet the boundary condition of no heat flow- 
ing out the crystal face. With our values of 5, £, a, P b , 
and Pi r , the calculated temperature T' agrees with our 
data to about one percent. 

It is interesting that the above treatment of BLIIRA 
does not distinguish between blue-light induced infrared 
absorption and infrared light induced blue-light absorp- 
tion. This model of BLIIRA only postulates that there 
is a two-color two-photon absorption cross-section. This 
treatment of two-color two-photon absorption in KNb03 
is perhaps simplistic. However, it serves the purpose in 
this paper of demonstrating that thermal self-locking of 
a Fabry-Perot cavity can be used to sensitively charac- 
terize the linear and nonlinear optical properties of an 
absorbing crystal. 
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VI. NEW MEASUREMENTS OF dn/dT 

From the calibration measurements that convert r\ in 
Figures | and | to a temperature rise, we can extract the 
temperature-dependent change in the index of refraction 
for KNDO3. When the crystal temperature changes, the 
optical pathlength of the bowtie cavity also changes be- 
cause both the crystal length and the index of refraction 
depend on temperature: 



dC _ 



dn x 

— + a T {n a 



no) 



(13) 



where C is the cavity optical pathlength, L c is the crystal 
length (10 mm), ax is the coefficient of linear expansion 
for KNbC>3, and n x and no are the indices of refraction 
for the crystal and air, respectively. 

In a relatively simple measurement, we can determine 
dn/dT for the crystal across a wide range of tempera- 
tures. We fix the laser frequency and the cavity mir- 
ror positions. We attenuate the laser intensity to 10/iW 
to avoid thermal self-locking. The crystal temperature 
can be independently controlled because it is mounted 
on a Peltier device. We measure the cavity transmission 
while slowly sweeping the crystal temperature over sev- 
eral degrees. At certain temperatures, the cavity optical 
pathlength meets the resonance condition, and we mea- 
sure a peak in the cavity transmission. The temperature 
between transmission peaks in these measurements corre- 
sponds to changing the optical pathlength by one optical 
wavelength, 846 nm. We repeat these measurements over 
a wide temperature range to determine dC/dT from -12 
°C to about 50 °C. 

Using equation |l3|, we convert dC/dT to dn/dT. With 
n c = 2.13, n b = 2.28, a T = 5 x 10~ 6 °CT 1 , and 
dC/dT = X/ AT, we determine dn/dT for several tem- 
peratures. Our values are shown in Figure |^ for light 
polarized along the c-axis and 6-axis, as labeled. Also 
shown in the figure is the data from Ghosh jhlj . The line 
from Ghosh is a fit to his measurements of the index of 
refraction between and 140 °C. We extrapolate that fit 
to lower temperatures to compare with our data. 

For light polarized along the c-axis, our measurements 
agree with those of Ghosh over our entire temperature 
range. For light polarized along the 6-axis, in the range 
5 to 50 °C, which is roughly the range for which our 
data overlaps Ghosh's data, our measurements also agree. 
However, below ~ 5°C, our data departs slightly from 
Ghosh's. 

To check our data, we made an additional measure- 
ment. When the polarization of the light incident on the 
cavity is rotated slightly relative to the c-axis, the trans- 
mission lineshape has two features, one for each polariza- 
tion. Because the index of refraction is different for the 
two polarizations (6 and c) , the optical pathlength is also 
different, and each polarization is transmitted at a dif- 
ferent cavity length. When we change the temperature 
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FIG. 6: The change in the index of refraction with temper- 
ature for light polarized along the b and c-axes. The points 
are from this work. The solid line is from Ghosh jyj. The 
dashed line is a second order polynomial fit to the data. 



of the crystal, the two peaks in the transmission line- 
shape shift by different amounts in different directions, 
as shown in Figure (7). We measured these shifts for a set 
of small temperature changes near -11.5 °C and deter- 
mined the relative changes in the indices of refraction, 
(dn c /dT) I (drib/dT), These measurements confirm our 
low temperature data in Figure [| Following the treat- 
ment of Ghosh, we fit our dn/dT data for the 6-axis to a 
second-order polynomial, and find 



drib a ^ 

= 1.55809 x 10" 9 T 2 

dl 

+ 4.06912 x 10~ 8 T - 33.76305 x 10" 



(14) 



where the units of dn/dT are °C , and T is measured 
in Celsius. 



VII. CONCLUSION 

We have demonstrated how thermal self-locking of a 
Fabry-Perot cavity can be used to determine optical 
properties of nonlinear crystals. We have determined 
new values for the one- and two-photon absorption co- 
efficients in KNbOs at 846 nm. We have also determined 
new values for the one photon (linear) absorption coef- 
ficient at 423 nm, and a BLIIRA absorption coefficient 
under specific circumstances. Finally, we present new 
measurements of the temperature-dependent change in 
the indices of refraction for the 6 and c axes. 

These measurements all derive from the thermal re- 
sponse of an absorbing Fabry-Perot cavity to light circu- 
lating in the cavity. The methods are general and can 
be widely applied to other systems, particularly those in 
which the absorption coefficients are small compared to 
scattering and reflection coefficients. 
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FIG. 7: Transmitted lineshape as a function of cavity length 
for two temperatures. The upper trace has been displaced for 
clarity. The polarization of the light incident on the cavity is 
rotated a few degrees relative to the c-axis. The larger feature 
is polarized along the c-axis, and the smaller feature is polar- 
ized along the &-axis, as labeled. This shows the magnitude 
and direction of the change in the cavity resonance length 
with temperature. 
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